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Abstract 

Presenting theoretical arguments and numerical results we demonstrate long-range intrachain 
correlations in concentrated solutions and melts of long flexible polymers which cause a systematic 
swelling of short chain segments. They can be traced back to the incompressibility of the melt 
leading to an effective repulsion u(s) ~ s/pR 3 (s) as Ce/i/s when connecting two segments together 
where s denotes the curvilinear length of a segment, R(s) its typical size, c e ~ V/^e ^ e "swelling 
coefficient", b e the effective bond length and p the monomer density. The relative deviation of 
the segmental size distribution from the ideal Gaussian chain behavior is found to be proportional 
to u(s). The analysis of different moments of this distribution allows for a precise determination 
of the effective bond length b e and the swelling coefficient c e of asymptotically long chains. At 
striking variance to the short-range decay suggested by Flory's ideality hypothesis the bond-bond 
correlation function of two bonds separated by s monomers along the chain is found to decay 
algebraically as 1/s 3 / 2 . Effects of finite chain length are considered briefly. 
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I. FLORY'S IDEALITY HYPOTHESIS REVISITED 



A cornerstone of polymer physics. Polymer melts are dense disordered systems con- 
sisting of macromolecular chains [1]. Theories that predict properties of chains in a melt 
or concentrated solutions generally start from the "Flory ideality hypothesis" formulated 
already in the 1940s by Flory 

aaa. 

This cornerstone of polymer physics states that 
chain conformations correspond to "ideal" random walks on length scales much larger than 
the _ dialer fl 1 fl Q. The _.v accepted justification of thi 8 m ean-field 
result is that mtrachain and interchain excluded volume forces compensate each other if 
many chains strongly overlap which is the case for three-dimensional melts |5|]. Since these 
systems are essentially incompressible, density fluctuations are known to be small. Hence, 
all correlations are supposed to be short-ranged as has been systematically discussed first 
by Edwards who developed the essential statistical mechanical tools , llO| also used 

in this paper. 

One immediate consequence of Flory's hypothesis is that the mean-squared size of chain 
segments of curvilinear length s = m — n (with 1 < n < m < N) should scale as -R 2 (s) = 
(r 2 ) = b 2 s if the two monomers n and m on the same chain are sufficiently separated along 
the chain backbone, and local correlations may be neglected (1 <C s). For the total chain 
(s = N-l > 1) this implies obviously that R 2 C {N - 1) = b 2 c (N - 1) « b 2 c N. Here, N denotes 
the number of monomers per chain, r the end-to-end vector of the segment, r = ||r|| its 
length and b e the "effective bond length" of asymptotically long chains [6|. (See Fig. [I] for 
an illustration of some notations used in this paper.) For the 2p-th moment (p = 0, 1, 2, . . .) 
of the segmental size distribution G(r, s) in three dimensions one may write more generally 
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which is, obviously, consistent with a Gaussian segmental size distribution 
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Both equations are expected to hold as long as the moment is not too high for a given 
segment length and the finite-extensibility of the polymer strand remains irrelevant [6]. 
Deviations caused by the segmental correlation hole effect. Recently, Flory's hypothe- 



sis has been challenged both theoretically [ll, [12|, [13|, \l4 , 
dimensional solutions 16, 0, Q, Q, |20j and ultrathin films 



14 . 05 1 and numerically for three- 
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22j | . These studies suggest 



that intra- and interchain excluded volume forces do not fully compensate each other on 
intermediate length scales, leading to long-range intrachain correlations. The general phys- 
ical idea behind these correlations is related to the "segmental correlation hole" of a typical 

n 

chain segment [19J. As sketched in Fig. [21 this induces an effective repulsive interaction when 
bringing two segments together, and swells (to some extent) the chains causing, hence, a 



systematic violation o 



briefly elsewhere 
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\. (CD . Elaborating and clarifying various points already presented 



20], we focus here on melts of long and flexible polymers. Us- 
ing two well-studied coarse-grained polymer models [23] various intrachain properties are 
investigated numerically as functions of s and compared with predictions from first-order 
perturbation theory. (For a discussion of intrachain correlations in reciprocal space see 

Refs. \m-) 

Central results tested in this study. The key claim verified here concerns the deviation 
5G(r,s) = G(r,s) — G (r,s) of the segmental size distribution G(r,s) from Gaussianity, 
Eq. (TJJ), for asymptotically long chains (N — > oo) in the "scale-free regime" (1 < s < N). 
We show that the relative deviation divided by c e / y/s scales as a function f(n) of n = r/b e \fs: 

8G(r,s)/G Q (r,s) _ ( 2 ^ 9 , 
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As we shall see, this scaling holds indeed for sufficiently large segment size r and curvilinear 
length s. The indicated "swelling coefficient" c e has been predicted analytically, 
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(p being the monomer number density), where we shall argue that the bond length of the 
Gaussian reference chain of the perturbation calculation must be renormalized to the effective 
bond length b e . Accepting Eq. (j3J) the swelling of the segment size is readily obtained by 
computing (r 2p ) = An J dr r 2+2p G(r, s). For the 2p-th moment this yields 

Kp{s) -2(2 P+ 1)\ yfs { ) 

For instance, for the second moment (p = 1) this reduces to Ki(s) = 1 — R e (s) 2 /b^s = 
c\l\fs ~ c e /\fs. We have replaced in Eq. (jSJ) the theoretically expected swelling coefficient 
c e by empirically determined coefficients c p . It will be shown, however, that c p /c e is close 
to unity for all moments. Effectively, this reduces Eq. (JSJ) to an efficient one-parameter 
extrapolation formula for the effective bond length b e of asymptotically long chains albeit 



empirical and theoretical swelling coefficients may slightly differ. While we show how b c may 
be fitted, no attempt is made to predict it from the operational model parameters and other 
measured properties such as the microscopic structure or the bulk compression modulus 



id- 



Outline. We begin our discussion by sketching the central theoretical ideas in Sec. [TT1 
There we will give a simple scaling argument and outline very briefly some elements of the 
standard perturbation calculations we have performed to derive them (Sec. IIIBj) . Details of 



the analytical treatment are relegated to Appendix |A] The numerical models and algorithms 
allowing the computation of dense melts containing the large chain lengths needed for a clear- 
cut test are presented in Sec. IIHI Our computational results are given in Sec. [TV] While 
focusing on long chains in dense melts, we explain also briefly effects of finite chain size. The 
general background of this work and possible consequences for other problems of polymer 
science are discussed in the final Sec. [V] 

II. PHYSICAL IDEA AND SKETCH OF THE PERTURBATION CALCULATION 
A. Scaling arguments 

Incompressibility and correlation of composition fluctuations. Polymer melts are essen- 
tially incompressible on length scales large compared to the monomer diameter, and the 
density p of all monomers does not fluctuate. On the other hand, composition fluctuations 
of labeled chains or subchains may certainly occur, however, subject to the total density 
constraint. Composition fluctuations are therefore coupled and segments feel an entropic 



penalty when their distance becomes comparable to their size [121 . Il9| . As sketched in 
Fig. E(a), we consider two independent test chains of length s in a melt of very long chains 
(N — > oo). If s is sufficiently large, their typical size, R(s) ~ b e y/s, is set by the effective 
bond length b e of the surrounding melt (taking apart finite chain-size effects). The test 
chains interact with each other directly and through the density fluctuations of the sur- 
rounding melt. The scaling of their effective interaction may be obtained from the potential 
of mean force U(r, s) = — ln(g(r, s)/g(oo, s)) where g(r) is the probability to find the second 
chain at a distance r assuming the first segment at the origin (r = 0). Since the correlation 
hole is shallow for large s, expansion leads to U(r, s) ~ 1 — g(r, s)/g(oo, s) ~ c(r, s)/p with 
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c(r, s) being the density distribution of a test chain around its center of mass. This distri- 
bution scales as c(r « 0, s) ~ s/R(s) d close to the center of mass (d being the dimension of 
space) and decays rapidly at distances of order R(s) paj]. Hence, the interaction strength at 
r/R(s) < 1 is set by u(s) = U(0, s) = c(0, s)/p « s/pR(s) d ~ s 1 "^ 2 Q, Q. Interestingly, 
■u(s) does not depend explicitly on the bulk compression modulus v. It is dimensionless and 
independent of the definition of the monomer unit, i.e. it does not change if A monomers 
are regrouped to form an effective monomer (p — > pj A, s —>■ s/X) while keeping the segment 
size R fixed. 

Connectivity and swelling. To connect both test chains to form a chain of length 2s 
the effective energy u(s) has to be paid and this repulsion will push the half-segments 
apart. We consider next a segment of length s in the middle of a very long chain. All 
interactions between the test segment and the rest of the chain are first switched off but we 
keep all other interactions, especially within the segment and between the segment monomers 
and monomers of surrounding chains. The typical size R{s) of the test segment remains 
essentially unchanged from the size of an independent chain of same strand length. If we 
now switch on the interactions between the segment and monomers on adjacent segments 
of same length s, this corresponds to an effective interaction of order u(s) as before. (The 
effect of switching on the interaction to all other monomers of the chain is inessential at 
scaling level, since these other monomers are more distant.) Since this repels the respective 
segments from each other, the corresponding subchain is swollen compared to a Gaussian 
chain of non-interacting segments. It is this effect we want to characterize. 

Perturbation approach in three dimensions. In the following we will exclusively consider 
chain segments s which are much larger than the number of monomers g = 1/vp contained 
in a blob [5], i.e. we will look on a scale where incompressibility matters. (The number g 
is also sometimes called "dimensionless compressibility" [ijj].) Interestingly, when taken at 
s = g the interaction strength takes the value 

with G z being the standard Ginzburg parameter used for the perturbation calculation 
of strongly interacting polymers |6j. Hence, the segmental correlation hole potential 
u(s) ~ G z (g /s) d / 2_1 G z for d > 2 and s ^> g. Although for real polymer melts as 
for computational systems large values of G z ~ 1 may sometimes be found, u(s) ~ 



decreases rapidly with s in three dimensions, as illustrated in Fig. MJd), and standard per- 
turbation calculations can be successfully performed. 

As sketched in the next paragraph these calculations yield quantities K[u] which are 
defined such that they vanish (K[u = 0] = 0) if the perturbation potential u(s) is switched 
off and are then shown to scale, to leading order, linearly with u. For instance, for the 
quantity K p (s), defined in Eq. ([1]), characterizing the deviation of the chain segment size 
from Flory's hypothesis one thus expects the scaling 

K p [u(s)} « +ti(a) « (7) 

The +- sign indicated marks the fact that the prefactor has to be positive to be consis- 
tent with the expected swelling of the chains. Consequently, the typical segment size, 
R(s)/b c ^ ml — u(s), must approach the asymptotic limit for large s from below. For 
three dimensional solutions Eq. (JTj) implies that K p (s) should vanish rapidly as l/(pbly/s). 
(This is different in thin films where u(s) ~ G z decays only logarithmically pjj as may be 
seen from Eq. fl33|) given below.) Taking apart the prefactors — which require a full calcu- 
lation — this corresponds exactly to Eq. (JSj) with a swelling coefficient c e m c p ~ 1/pb^. in 
agreement with Eq. (j3J). Note also that the predicted deviations are inversely proportional 
to 6g, i.e. the more flexible the chains, the more pronounced the effect. Similar relations 
K[u] ~ u may also be formulated for other quantities and will be tested numerically in 
Sec. IIVI There, we will also check that the linear order is sufficient. 

B. Perturbation calculation 

Generalities. Before delving more into our computational results we summarize here how 
Eqs. (J3JISJ) and related relations have been obtained using standard one-loop perturbation 
calculation. The general task is to determine (A) ~ (*4.) (1 + (U) Q ) — (AU) for measurable 
quantities A such as the squared distance between two monomers n and m on the same chain, 
A = r^ m . Here, (...) denotes the average over the distribution function of the unperturbed 
ideal chain of bond length b and U = dk dlv(rki) the effective perturbation potential. 
We discuss first the general results in the scale free regime (1 < s C N), argue then that b 
should be renormalized to the effective bond length b e and sketch finally the calculation of 
finite chain-size effects. 
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The scale free regime. Following Edwards [61,3,0], the Gaussian (or "Random Phase" 
[5I) approximation of the pair interaction potential in real space is 



where v is a parameter which tunes the monomer interaction. (It is commonly associated 
with the bare excluded volume of the monomers 6j, but should more correctly be identified 
with the bulk modulus effectively measured for the system. See the discussion of Eq. (15) of 
Ref. ([3]).) The effective potential consists of a strongly repulsive part v5(r) of very short 
range, and a weak attractive part of range £ where the correlation length of the density 
fluctuations is given by £ 2 = b 2 g/12 with g = 1/pv. In Fourier space Eq. ([8]) is equivalent to 

q 2 

with q being the wave vector. This is sufficient for calculating the scale free regime corre- 
sponding to asymptotically long chains where chain end effects may be ignored. The different 
graphs one has to compute are indicated in Fig. [TJ For A = r 2 (with 1 < n < m C Af) this 
yields, e.g., 



(r) = b {1 + — — j(m-n) 



n 



«»|i-t=(t))- '■"» 
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In the second line we have used the definition of the swelling coefficient c e indicated in 
Eq. fll]) and have set 

« s K 1+i v£)=*( i+ ^ G -) 

with G z = ^Jvp/b z p and p — 1. (The prefactor p has been added for convenience.) The 
coefficient b e of the leading Gaussian term in Eq. flTUJ) - - entirely due to the graph Jj 
describing the interactions of monomers inside the segment between n and m — has been 

n 

predicted long ago by Edwards [6]. It describes how the effective bond length is increased 
from b to b c under the influence of a small excluded volume interaction. The second term in 
Eq. (TTDT) entails the l/-\/i-swelling which is investigated numerically in this paper. It does 
only depend on b and p but, more importantly, not on v — in agreement with the scaling 
of u(s) discussed in Sec. Ill Al The relative weights contributing to this term are indicated 



in Fig. [TJin units of — ^ Q / tt 3 vS, 2 / b 3 y/s . The diagrams J_ and I + are obviously identical in 
the scale free limit. Note that the interactions described by the strongest graph Jj align the 
bonds l n and l m while the others tend to reduce the effect. 

For higher moments of the segment size distribution G(r, s) it is convenient to calculate 
first the deviation of the Fourier-Laplace transformation of 5G(r, s) and to obtain the mo- 
ments from the coefficients of the expansion of this "generating function" in terms of the 
squared wave vector q 2 . As explained in detail in the Appendix |A] this yields more generally 

(r « _ (fr+py r d _ sew * ( b \"A (12) 



where we have used Eq. (iTTj) with general p. Obviously, Eq. (1121) is consistent with our 
previous finding Eq. ffTUl) for p = 1. The corresponding segmental size distribution is 

G(r ' s) = Gi) exp H£ 

3 \ 3/2 / 3r 2 



eX H-2^J^UJ /W (13) 

with n = r/hyfs and /(n) being the same function as indicated in Eq. ([3]). The leading 
Gaussian terms in Eqs. (TT2]) and (ITBl depend on the effective bond length 6 e , the second 
only on the Kuhn length b of the reference chain. When comparing these result with Eqs. (J3j) 
and (jSJ) proposed in the Introduction, one sees that both equations are essentially identical - 
taken apart, however, that they depend on b and b e . Note the conspicuous factor (b/b e ) 2p ~ 3 
in Eq. f|T2|) which would strongly reduce the empirical swelling coefficients c p = c e (b/b e ) 2p ~ 3 
for large p if b and b e were different. 

Interpretation of first-loop results in different contexts. The above perturbation results 
may be used directly to describe the effect of a weak excluded volume v on a reference system 
of perfectly ideal polymer melts with Kuhn segment length b where all interactions have been 
switched off (v — 0). It is expected to give a good estimation for the effective bond length 
b e only for a small Ginzburg parameter: G z < 1. For the dense melts we want to describe 
this does not hold (Sec. IIII[) and one cannot hope to find a good quantitative agreement 
with Eq. (ITT]) . Note also that large wave vectors contribute strongly to the leading Gaussian 
term. The effective bond length b e is, hence, strongly influenced by local and non-universal 
effects and is very difficult to predict in general. 
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Our much more modest goal is to predict the coefficient of the 1 / y^-perturbation and 
to express it in terms of a suitable variational reference Hamiltonian characterized by a 
conveniently chosen Kuhn segment b and the measured effective bond length b e (instead of 
Eq. (jlip ). Following Muthukumar and Edwards [loj], we argue that for dense melts b should 
be renormalized to b c to take into account higher order graphs. No strict mathematical 
proof can be given at present that the infinite number of possible graphs must add up in 
this manner. Our hypothesis relies on three observations: 

• The general scaling argument discussed in Sec. Ill Al states that we have only one 
relevant length scale in this problem, the typical segment size R(s) ~ b e ^fs itself. The 
incompressibility constraint cannot generate an additional scale. It is this size R{s) 
which sets the strength of the effective interaction which then in turn feeds back to 
the deviations of R(s) from Gaussianity. Having a bond length b in addition to the 
effective bond length b e associated with R(s) would imply incorrectly a second length 
scale b\fs varying independently with the bulk modulus v. (We will check explicitly 
below in Fig. [13] that there is only one length scale.) This implies b/b e = const v°. 

• Thus, since by construction b/b e = 1 for v — > 0, it follows that both lengths should be 
equal for all v. 

• We know from Eq. (IT2I) that the empirical coefficients c p = c e (b/b e ) 2p ~ 3 should depend 
strongly on the moment considered if the ratio b/b e is not close to unity. It will be 
shown below (Fig. [6]) that c p /c e ~ 1 for all p. This implies b ~ b e . 

Finite chain size effects. To describe properly finite chain size corrections Eq. (Q must 
be replaced by the general linear response formula 

-+ P F(q) (14) 



v(q) v 

with F(q) = Nf£,(x) being the form factor of the Gaussian reference chain given by Debye's 
function fo{x) = 2(e~ x — 1 + x)/x 2 with x = (qb) 2 N/6 [fj. This approximation allows 
in principle to compute, for instance, the (mean-squared) total chain end-to-end distance, 
A = (rjv — ri) 2 . One verifies readily (see |6[], Eq. (5. III. 9)) that the effect of the perturbation 
may be expressed as 

d 3 q q 2 b 4 /"^ , / q 2 b 2 s 



(A} (U} - (AU} = J %) J ds S ^ N _ s) exp 
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(15) 



We take now first the integral over s. In the remaining integral over q small q wave vectors 
contribute to the A/iV-swelling while large q renormalize the effective bond length of the 
dominant Gaussian behaviour linear in N (as discussed above). Since we wish to determine 
the non-Gaussian corrections, we may focus on small wave vectors q <C l/£. Since in this 
limit 1/v = pg <C pF(q), one can neglect in Eq. ( 1T41) the 1/v contribution to the inverse 
effective interaction potential. We thus continue the calculation using the much simpler 
v(N,x) = l/(Npf D (x)). This allows us to express the swelling as 

blN [uh 1 j 

To simplify the notation we have set here finally b = b c in agreement with the hypothesis 
discussed above. The numerical integral I(x u ) = J* u dx . . . over x is slowly convergent at 
infinity. As a consequence the estimate I{oo) = 1.59 may be too large for moderate chain 
lengths. In practice, convergence is not achieved for values x u (N) « (b/£/) 2 N corresponding 
to the screening length £. 

We remark finally that numerical integration can be avoided for various properties if the 
Pade approximation of the form factor, F(q) = N/(l + (qb) 2 N/12), is used. This allows 
analytical calculations by means of the simplified effective interaction potential 



q 



b 5 b 3 



f '<«) = 1^5 + iV^- < 17) 
This has been used for instance for the calculation of finite chain size effects for the bond- 



bond correlation function discussed in Sec. IIVGI below 6l| 



III. COMPUTATIONAL MODELS AND TECHNICAL DETAILS 
A. Bond fluctuation model 

A widely-used lattice Monte Carlo scheme for coarse-grained polymers. The body of 
our numerical data comes from the three dimensional bond fluctuation model (BFM) - 
a lattice Monte Carlo (MC) algorithm where each monomer occupies eight sites of a unit 



cell of a simple cubic lattice 



24 



261 ] . Our version of the BFM with 108 bond vectors 



corresponds to flexible athermal chain configurations 23j. All length scales are given in units 
of the lattice constant and time in units of Monte Carlo Steps (MCS). We use cubic periodic 
simulation boxes of linear size L = 256 containing n mon = pL 3 = 2 20 w 10 6 monomers. 
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This monomer number corresponds to a monomer number density p = 0.5/8 where half of 
the lattice sites are occupied (volume fraction 0.5). The large system sizes used allow us to 
suppress finite box-size effects for systems with large chains. Using a mix of local, slithering 



snake 



28, 



29| , and double-bridging 



13, 



23, 



30 



3JJ MC moves we were able to equilibrate 



dense systems with chain lengths up to iV = 8192. 

Equilibration and sampling of high-molecular BFM melts. Standard BFM implementa- 



tions 



26 



32 



use local MC jumps to the 6 closest lattice sites to prevent the crossing 
of chains and conserve therefore the chain topology. These "L06" moves lead to very large 
relaxation times, scaling at least as r c ~ iV 3 , as may be seen from Fig. [3] (stars). The relax- 
ation time r e = R^/6D S indicated in this figure has been estimated from the self-diffusion 
coefficient D s obtained from the mean-square displacements of all monomers in the free 
diffusion limit. (For the largest chain indicated for L06 dynamics only a lower bound for 
r e is given.) Instead of this more realistic but very slow dynamical scheme we make jump 
attempts to the 26 sites of the cube surrounding the current monomer position (called "L26" 
moves). This allows the chains to cross each other which dramatically speeds up the dynam- 
ics, especially for long chains (N > 512). If only local moves are considered, the dynamics 
is perfectly consistent with the Rouse model {f| . As shown in Fig. [31 we find r e w 530 iV 2 for 
L26 dynamics. This is, however, still prohibitive b y la rge for sampling configurations with 



62|. 



the longest chain length N we aim to characterize 

Slithering snake moves. In addition to the local moves one slithering snake move per 
chain is attempted on average per MCS corresponding to the displacement of N monomer 
along the chain backbone. Note that in our units two spatial displacement attempts per 
MCS are performed on average per monomer, one for a local move and one for a snake 
move. (In practice, it is computationally more efficient for large N to take off a monomer 
at one chain end and to paste it at the other leaving all other monomers unaltered. Before 
dynamical measurements are performed the original order of beads must then be restored.) 
Interestingly, a significantly larger slithering snake attempt frequency would not be useful 



since the relaxation time o: 
exponentially with mass 29 



slithering snakes without or only few local moves increases 



34| due to the correlated motion of snakes 35j. In order to 



obtain an efficient free snake diffusion (with a chain len gth independent curvilinear diffusion 



coefficient D C (N) ~ iV° and r c « N 2 /D C (N) ~ iV 2 [2J, [29|]) it is important to relax density 



fluctuations rapidly by local dynamical pathways. As shown in Fig. [3] (squares), we find 
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a much reduced relaxation time r e ~ 40N 2 which is, however, still unconveniently large 
for our longest chains. Note that most of the CPU time is still used by local moves. The 
computational load per MCS remains therefore essentially chain length independent. 

Advantages and pitfalls of double-bridging moves. Double-bridging (DB) moves are very 
useful for high densities and help us to extend the accessible molecular masses close to 
10 4 . As for slithering snake moves we use all 108 bond vectors to switch chain segments 
between two different chains. Only chain segments of equal length are swapped to conserve 
monodispersity. Topolocial constraints are again systematically and deliberately violated. 
Since more than one swap partner is possible for a selected first monomer, delicate detailed 
balance questions arise. This is particularly important for short chains and is discussed 



in detail in Ref. 



23]. Technically, the simplest solution to this problem is to refuse all 



moves with more than one swap partner (to be checked both for forward and back move). 
The configurations are screened with a frequency /db for possible DB moves where we 
scan in random order over the monomers. The frequency should not be too large to avoid 
(more or less) immediate back swaps and monomers should move at least out of the local 
monomer cage and over a couple of lattice sites. We use /db = 0.1 between DB updates for 
the configurations reported here. (The influence of /db cm the performance has not been 
explored systematically, but preliminary results suggest a slightly smaller DB frequency for 
future studies.) The diffusion times over the end-to-end distance for this case are indicated 
in Tab. B As shown in Fig. El we find empirically r c (N) « 13iV L62 . For N = 8192 this 
corresponds to 3 • 10 7 MCS. This allows us even for the largest chain lengths to observe 
monomer diffusion over several R e within the 10 8 MCS which are feasible on our XEON-PC 
processor cluster. 

The efficiency of DB moves is commonly characterized in terms of the relaxation time 
r ee of the end-to-end vector correlation function [30I . 31]. For normal chain dynamics this 



would indeed characterize the longest relaxation time of the system, i.e. r e ~ r ee . For the 
double-bridging this is, however, not sufficient since density fluctuations do not couple to the 
bridging moves and can not be relaxed. We find therefore that configurations equilibrate 
on time scales given by r c rather than by r ce <C r c . This may be verified, for instance, 
from the time needed for the distribution R e (s) (and especially its spatial components) to 



equilibrate. The criterion given in the literature 



30] is clearly not satisfactory and may lead 



to insufficiently equilibrated configurations. In summary, equilibration with DB moves still 
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requires monomer diffusion over the typical chain size, however at a much reduced price. 

Some properties of our configurations. The Tables [J and [III summarize some system 
properties obtained for our reference density p = 0.5/8. Averages are performed over all 
chains and 1000 configurations. These configurations may be considered to be independent 
for iV < 4096. Only a few independent configurations exist for the largest chain length 
N = 8192 which has to be considered with some care. Taking apart this system, chains are 
always much smaller than the box size. For asymptotically long chains, we obtain an average 
bond length (|1|) rs 2.604, a root-mean-squared bond length I = (l 2 ) 1 ^ 2 « 2.635 and an 
effective bond length b e ~ 3.244 — as we will determine below in Sec. IIV Al This corresponds 
;o a ratio = h\jl 2 « 1.52 and, hence, to a persistence length l p = I (Coo + l)/2 ~ 3.32 



23j. Especially, we find from the zero wave vector limit of the total structure factor S(q) 



a low (dimensionless) compressibility g = S(q — ► 0)/p ~ 0.246 which compares well with 
real experimental melts. From the measured bulk compression modulus v = l/g(p)p and 
the effective bond length b e one may estimate a Ginzburg parameter G z = ^fvpjo\p ~ 0.96. 
Following Ref. [l^] the interaction parameter v is supposed here to be given by the full 
inverse compressibility and not just by the second virial coefficient. 



B. Bead spring model 



Hamiltonian. Additionally, molecular dynamics simulations of a bead-spring model 
(BSM) 36j were performed to dispel concerns that our results are influenced by the un- 
derlying lattice structure of the BFM. The model is derived from a coarse-grained model 
for polyvinylalcohol which has been employed to study polymer crystallization [3^]. It is 
characterized by two potentials: a non-bonded potential of Lennard- Jones (LJ) type and a 



harmonic bond potential. While the often employed Kremer-Grest model |3§J uses a 12 — 6 
LJ potential to describe the non-bonded interactions U n h(r), our non-bonded potential has 
a softer repulsive part. It is given by 



= 1.511 p"_ 
which is truncated and shifted at the minimum at r mm w 



(18) 



1.15. Note that all length scales 



are given in units of o"o and we use LJ units 39j for all BSM data (mass m = 1, Boltzmann 
constant k B = 1). The parameters of the bond potential, U\,(r) = 1120(r — Z b ) 2 , are adjusted 
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so that the average bond length Up = 0.84) « £ b = 0.97 is approximately the same as in the 
standard Kremer-Grest model [38(] . The average bond length and the root-mean-squared 
bond length are almost identical for the BSM due to the very stiff bond potential. Since 
r minA ~ 1.16 bonded monomers penetrate each other significantly. 

Equilibration and sampling. We perform standard molecular dynamics simulations in the 
canonical ensemble with a Langevin thermostat (friction constant T = 0.5) at temperature 



T — 1. The equations of motion are integrated by the velocity- Ver let algorithm 39]. To 
improve the statistics for large chain length, we have implemented additional double-bridging 
moves. Since only few of these MC moves are accepted per unit time, this does affect neither 
the stability nor the accuracy of the molecular dynamics sweeps. 

Some properties obtained. For clarity, we show only data for chain length iV = 1024 
and number density p = 0.84, the typical melt density of the Kremer-Grest model [38]. 
For the reported data we use periodic simulation boxes of linear size L m 62 containing 
n mon = 196608 monomers, but we have also sampled different boxes sizes (up to L — 77.5) 
to check for finite box-size effects. For the reference density a dimensionless compressibility 
g ~ 0.08 is found which is about three times smaller than for our BFM melt. For the 
effective bond length we obtain b e « 1.34, i.e. BSM chains (C^ ~ 1.91, l p ~ 1.41) are 
slightly stiffer than the corresponding BFM polymers. Fortunately, the product pb^ ~ 2 is 
roughly similar in both models and one expects from Eq. (T4|) a similar swelling for large 
s. Note finally that the Ginzburg parameter G z ~ 1.8 is much larger than for the BFM 
systems. As we have emphasized in Sec. [TTJ, this should, however, not influence the validity 
of the perturbation prediction of the expected 1 / -^/i-swelling of the chains when expressed 
in terms of the measured effective bond length. 



IV. NUMERICAL RESULTS 



As illustrated in Fig. [TJ a chain segment of curvilinear length s > is identified by two 
monomers n and m = n + s on the same chain. We compute here various moments of 
chain segment properties where we ensemble-average over all chains and all start points n. 
The statistical accuracy must therefore always decrease for large s. We concentrate first on 
the second moment (p = 1) of the segmental size distribution. Higher moments and the 
segmental size distribution are discussed in Sec. IIVEI 
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A. The swelling of chain segments 

Scale free regime for 1 < s < ^. The mean-squared segment size R^{s) = (r 2 ) is 
presented in the Figs. HJO and The first plot shows clearly that chain segments are swollen, 
i.e. R^(s)/s increases systematically and this up to very large curvilinear distances s. Only 
BFM data are shown for clarity. A similar plots exists for the BSM data. In agreement 
with Eq. (j3j) for p — 1, the asymptotic Gaussian behavior (dashed line) is approached 
from below and the deviation decays as u(s) oc 1/y/s (bold line). The bold line indicated 
corresponds to b c = 3.244 and c\ ~ c e ~ 0.41 which fits nicely the data over several decades 
in s — provided that chain end effects can be neglected (s <C A). Note that a systematic 
underestimation of the true effective bond length would be obtained by taking simply the 
largest Rl(s)/s ~ 3.23 2 value available, say, for monodisperse chains of length A = 2048. 

Finite chain-size effects. Interestingly, R%(s)/s does not approach the asymptotic limit 
monotonicly. Especially for short chains one finds a non-monotonic behavior for s —> N. 
This means that the total chain end-to-end distance R e (s — N — 1) must show even more 
pronounced deviations from the asymptotic limit. This is confirmed by the dashed line 
representing the b^(N) = R^(N — 1)/(A — 1) data points given in Tab. [H We emphasize 
that the non-monotonicity of i? 2 (s)/s becomes weaker with increasing A and that, as one 
expects, the inner distances, as well as the total chain size, are characterized by the same 
effective bond length b e for large s or A. The non-monotonic behavior may be qualitatively 
understood by the reduced self-interactions at chain ends which lessens the swelling on these 
scales. These finite- A corrections have been calculated analytically using the full Debye 
function for the effective interaction potential v(q), Eq. ( |T4T) . The prediction for the total 
chain end-to-end vector given in Eq. ffT6l) is indicated in Fig. [4] (dash-dotted line) where we 
have replaced the weakly A-dependent integral I(x u ) by its upper bound value for infinite 
chains 

, Re(N-l) 1.59c e 

b 2 c (N-l) y/N=T' { ] 

We have changed here the chain length A in the analytical formula (obtained for large chains 
where A pa A — 1 ) to the curvilinear length A — 1 . This is physically reasonable and allows 
to take better into account the behavior of small chains. Note that Eq. f|T9l) is similar to 
Eq. ((SD — apart from a slightly larger prefactor explaining the observed stronger deviations. 
Theory compares well with the measured data for large A. It does less so for smaller A, 
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as expected, where the chain length dependence of the numerical integral I(x u (N)) < 1.59 
must become visible. This explains why the data points are above the dash-dotted line. 
Note also that additional non-universal finite- N effects not accounted for by the theory are 
likely for small N. In contrast to this, R e (s) is well described by the theory even for rather 
small s provided that N is large and chain end effects can be neglected. In summary, it is 
clear that one should use the segment size R e {s) rather than the total chain size to obtain 
in a computational study a reliable fit of the effective bond length b e . 

Extrapolation of the effective bond length of asymptotically long chains. The represen- 
tation chosen in Fig. 0] is not the most convenient one for an accurate determination of b c 
and c\. How precise coefficients may be obtained according to Eq. (pxj) is addressed in the 
Figs. and The fitting of the effective bond length b e and its accuracy is illustrated in 
Fig. [5] for BFM chains of length N = 2048. This may be first done approximately in linear 
coordinates by plotting R%(s)/s as a function of l/y/s (not shown). Since data for large s 
are less visible in this representation, we recommend for the fine-tuning of b e to switch then 
to logarithmic coordinates with a vertical axis y — 1 — R^(s)/b^s for different trial values 
of b e . The correct value of b e is found by adjusting the vertical axis y such that the data 
extrapolates linearly as a function of 1/ \fs to zero for large s. We assume for the fine-tuning 
that higher order perturbation corrections may be neglected, i.e. we take Eq. (jSJ) literally. 
(We show below that higher order corrections must indeed be very small.) The plot shows 
that this method is very sensitive, yielding a best value that agrees with the theory over 
more than one order of magnitude without curvature. As expected, it is not possible to 
rationalize the numerically obtained values b e ~ 3.244 for the BFM and b e ~ 1.34 for the 
BSM using Eq. (Ilip . According to Eq. (T4J) these fit values imply the theoretical swelling 
coefficients c e = 0.41 for the BFM and c e = 0.44 for the BSM. 

Empirical swelling coefficients. As a next step the horizontal axis is rescaled such that 
all data sets collapse on the bisection line, i.e. using Eq. ([5]) we fit for the empirical swelling 
coefficient c\ and compare it to the predicted value c e . This rescaling of the axes allows to 
compare both models in Fig. El For clarity the BSM data have been shifted upwards. For 
the BFM we find ci/c e ~ 1.0, as expected, while our BSM simulations yield a slightly more 
pronounced swelling with ci/c e ~ 1.2. 

Segmental radius of gyration. Also indicated in Fig. [6] is the segmental radius of gyration 
R g (s) (filled circles) computed as usual jf| as the variance of the positions of the segment 
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monomers around their center of mass. Being the sum over all s + 1 monomers, it has a 
much better statistics compared to R e (s). The scaling used can be understood by expressing 
the radius of gyration Rg(s) = , ^ 2 (Y^k=n YM=n r li) m terms of displacement vectors r k i 
. Using Eq. (jHJ) and integrating twice this yields 



6R 2 g {s) _ 8 a 



(20) 



b 2 e (s + l) 5^s + l 

Plotting the l.h.s. of this relation against the r.h.s. we obtain a perfect data collapse on the 
bisection line where we have used the same parameters b e and C\ as for the mean-squared 
segment size. This is an important cross-check which we strongly recommend. Different 
values indicate insufficient sample equilibration. 



B. Chain connectivity and recursion relation 

As was emphasized in Sec. Ill Al the observed swelling is due to an entropic repulsion 
between chain segments induced by the incompressibility of the melt. To stress the role of 
chain connectivity we repeat the general scaling argument given above in a form originally 
proposed by Semenov and Johner for ultrathin films [3]. As shown in Fig. [7] we test the 
relation 

A( ) - (X-VX)Rl(s) U[S) ~ pRM d ( } 

with K\(s) being a direct measure of the non-Gaussianity (A being a positive number) 

comparing the size of a segment of length As with the size of A segments of length s joined 
together. (The prefactor 1/(A — VT) in the definition of K\(s) has been introduced for 
convenience.) Equivalently, this can be read as a measure for the swelling of a chain where 
initially the interaction energy u between the segments has been switched off. K\ is a 
functional of u(s) with K\\u = 0] = 0. The analytic expansion of the functional must be 
dominated by the linear term (as indicated by ~ in the above relation) simply because u is 
very small. Altogether, Eq. (1211) yields a recursion relation relating R e (Xs) with R e (s) for 
any A provided 1 s < Xs N. It can be solved, leading (in lowest order) to Eq. (J5j) with 
p = 1. This may be seen from the ansatz Rl{s) = bls(l — c e / + . . .) which readily yields 
uj = 3/2 and c e ~ l/p&e- 

Eq. ( !2T|) has been validated directly in Fig. [7] for A = 2 (corresponding to two segments 
of length s joined together) for the BFM and the BSM as indicated. In addition, for BFM 
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chains of length N = 2048 several values of A have been given. As suggested by Eq. (j3j), 
we have plotted K\(s) as a function of (ci/c c ) u(s) with u(s) = ^/2A/rc 3 s/ pRl(s) ~ c e /^/s. 
The prefactor of w(s) allows a convenient comparison with Fig. [5] Note the perfect data 
collapse for all data sets. More importantly, the predicted linearity is well confirmed for 
large segments (1 C s) and this without any tunable parameter for the vertical axis, as was 
needed in the previous Figs. [5] and El 



C. Intrachain bond-bond correlations 



Expectation from Flory's hypothesis. An even more striking violation of Flory's ideality 
hypothesis may be obtained by computing the bond-bond correlation function, defined by 
the first Legendre polynomial P(s) = (l OT = n +s ■ K) A 2 where the average is performed, as 
before, over all possible pairs of monomers (n, m = n + s) 63|]. Here, L. = r,; + i — denotes 
the bond vector between two adjacent monomers % and % + 1 and I 2 = (l 2 l ) n the mean- 
squared bond length. The bond-bond correlation function is generally believed to decrease 
exponentially [4]. This belief is based on the few simple single chain models which have 
been solved rigorously 4, 4^] and on the assumption that all long range interactions are 
negligible on distances larger than the screening length £. Hence, only correlations along 
the backbone of the chains are expected to matter and it is then straightforward to work 
out that an exponential cut-off is inevitable due to the multiplicative loss of any information 
transferred recursively along the chain 4]. 

Asymptotic behavior in the melt. That this reasoning must be incorrect follows imme- 
diately from the relation 

p < s > = WJ? R ^ (22) 

expressing the bond-bond correlation function as the curvature of the second moment of 
the segment size distribution. It is obtained from the identity (l n ■ l m ) = (<9 n r n ■ d m r m ) = 
—d n d m (r^ m ) /2. (Note that the velocity correlation function is similarly related to the 
second derivative of the mean-square displacement with respect to time [4l|.) Hence, P(s) 
allows us to probe directly the non-Gaussian corrections without any ideal contribution. 
This relation together with Eq. (jSJ) suggests an algebraical decay P(s) = c-p/s^ with 

cu = 3/2 , c P = ci (b c /l) 2 /8 w (23) 
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of the bond-bond correlation function for dense solutions and melts, rather than the ex- 
ponential cut-off expected from Flory's hypothesis. This prediction (bold line) is perfectly 
confirmed by the larger chains (N > 256) indicated in Fig. [HJ In principle, the swelling 
coefficient, C\ ~ cp, may also be obtained from the power law amplitude of the bond-bond 
correlation function, however, to lesser accuracy than by the previous method (Fig. One 
reason is that P(s) decays very rapidly and does not allow a precise fit beyond s ~ 10 2 . The 
values of cp obtained from c\ are indicated in Tab. [Til Data from the BSM have also been 
included in the figure to demonstrate the universality of the result. The vertical axis has 
been rescaled with cp which allows to collapse the data of both models. 

Finite chain-size corrections. As can be seen for N = 16, exponentials are compatible 
with the data of short chains. This might explain how the power law scaling has been 
overlooked in previous numerical studies, since good statistics for large chains (N > 1000) 
has only become available recently. However, it is clearly shown that P(s) approaches 
systematically the scale free asymptote with increasing N. The departure from this limit is 
fully accounted for by the theory if chain end effects are carefully considered (dashed lines). 
Generalizing Eq. ( f23l) and using the Pade approximation, Eq. ( 1TT1) . perturbation theory 
yields 

where we have set x = v/s/iV. For x C 1 this is consistent with Eq. (1231) . In the limit of 
large s — > N, the correlation functions vanish rigorously as P(s) oc (1 — x) 2 . Considering 
that non-universal features cannot be neglected for short chain properties and that the 
theory does not allow for any free fitting parameter, the agreement found in Fig. [S]is rather 
satisfactory. 



D. Higher moments and associated coefficients 

Effective bond length and empirical swelling coefficients. The preceding discussion fo- 
cused on the second moment of the segmental size distribution G(r, s). We have also com- 
puted for both models higher moments (r 2p ) with p < 5. If traced in log-linear coordinates 
as y — (6 p p\ (r 2p ) /{2p + l)\s p ) l l p vs. x = s higher moments approach b\ from below - 
just as the second moment presented in Fig. HI The deviations from ideality are now more 
pronounced and increase with p (not shown). The moments are compared in Fig. [6] with 
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Eq. (j3j) where they are rescaled as y — K p (s) as defined in Eq. ([T]) and plotted as functions 



°f x = 2(2p+i)! ^ ne prediction is indicated by bold lines. It is important that the 
effective bond length b e is obtained from the analysis of all functions K p (s) as illustrated in 
Fig. El Otherwise we would regard equilibration and statistics as insufficient. 

The empirical swelling coefficients c p are obtained, as above in Sec. IIV Al by shifting the 
data horizontally. A good agreement with the expected c p /c e « 1 is found for both models 
and all moments as may be seen from Tab. [Til This confirms the renormalization of the 
Kuhn segment b — > b e of the Gaussian reference chain in agreement with our discussion in 
Sec. Ill Bt Otherwise we would have measured empirical coefficients decreasing strongly as 
c p / c e ~ (b/b e ) 2p ~ 3 with p. Since the effective bond length of non-interacting chains are known 
for the BFM (b ~ 2.688) and the BSM (6 « 0.97), one can simply check, say for p = 5, that 
the non-renormalized values would correspond to the ratios c^/c e ~ (2.688/3.244) 7 « 0.3 for 
the BFM and c 5 /c e « (0.97/1. 34) 7 « 0.1 for the BSM. This is clearly not consistent with 
our data. 

It should be emphasized that both coefficients b e and c p are more difficult to determine 
for large p, since the linear regime for x <C 1 in the representation chosen in Fig. [6] becomes 
reduced. For large x ^> 1 one finds that y(x) — > 1, i.e. (r 2p ) jb 2p s v — > 0. This trivial 
departure from both Gaussianity and the l/y^-deviations we try to describe, is due to the 
finite extensibility of chain segments of length s which becomes more marked for larger 
moments probing larger segment sizes. The data collapse for both x-regimes is remarkable, 
however. Incidentally, it should be noted that for the BSM the empirical swelling coefficients 
are slightly larger than expected. At present we do not have a satisfactory explanation for 
this altogether minor effect, but it might be attributed to the fact that neighbouring BSM 
beads along the chain strongly interpenetrate — an effect not considered by the theory. 

Non- Gaussian parameter a p . The failure of Flory's hypothesis can also be demonstrated 
by means of the standard non-Gaussian parameter 

gp | / 2p \ 

„ / „\ i \j jj- \ 1 -nm/ fnK\ 

a ' {s> = 1 - WTTy. Taj (25) 

comparing the 2p-th moment with the second moment (p = 1). In contrast to the closely 
related parameter K p (s) this has the advantage that here two measured properties are com- 
pared without any tuneable parameter, such as b e , which has to be fitted first. Fig. presents 
a p (s) vs. c e / f° r f ne three moments with p = 2,3, 4. For each p we find perfect data col- 
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same 



lapse for all chain lengths and both models and confirm the linear relationship a p {s) ~ u(s) 
expected. The lines indicate the theoretical prediction 

3(2^)' \* 

2 (2p+l)\ V 1 ; 

which can be derived from Eq. ([5]) by expanding the second moment in the denominator. An 
alternative derivation based on the coefficients of the expansion of the generating function 
G(q, s) in q 2 is indicated by Eq. ( 1A2I) in the Appendix. Having confirmed above that c p /c e ~ 
1, we assume in Eq. ( 1261) that c p = c e to simplify the notation. The prefactors 6/5, 111/35 
and 604/105 for p = 2, 3 and 4 respectively are nicely confirmed. They increase strongly 
with p, i.e. the non-Gaussianity becomes more pronounced for larger moments as already 
mentioned. Note also the curvature of the data at small s due to the finite extensibility 
of the segments which becomes more marked for higher moments. If one plots az p (s) as a 
function of the r.h.s. of Eq. ( 1261) all data points for all moments and even for too small s 
collapse on one master curve (not shown) — just as we have seen before in Fig. (J6j). 

Correlations of different directions. A similar correlation function is presented in Fig. [10] 
which measures the non-Gaussian correlations of different spatial directions. It is defined by 

W.) - i - (27) 

for the two spatial components x and y of the vector r as illustrated by the sketch given 
at the bottom of Fig. [TUJ Symmetry allows to average over the three pairs of directions 
(x, y), (x,z) and (x,z). Following the general scaling argument given in Sec. [Ill we expect 
K xy (s) ~ u(s) ~ c e /\/i which is confirmed by the perturbation result 

K xy {s) = \^ = K 2 (s). (28) 

This is nicely confirmed by the linear relationship found (bold line) on which all data from 
both simulation models collapse perfectly. The different directions of chain segments are 
therefore coupled. As explained in the Appendix (Eq. (1A3j) ). K xy (s) and a 2 (s) must be 
identical if the Fourier transformed segmental size distribution G(q, s) can be expanded 
in terms of q 2 and this irrespective of the values the expansion coefficients take. Fig. [TUl 
confirms, hence, that our computational systems are perfectly isotropic and tests the validity 
of the general analytical expansion. 

The correlation function K xy is of particular interest since the zero-shear viscosity should 
be proportional to (o"^) ~ (x 2 y 2 ) = (x 2 ) (y 2 ) (1 — K xy (s)). We assume here following 
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Edwards [6] that only intrachain stresses contribute to the shear stress a xy . Hence, our 
results suggest that the classical calculations Q — assuming incorrectly K xy = — should 
be revisited. 



E. The segmental size distribution 

We turn finally to the segmental size distribution G(r, s) itself which is presented in 
Figs. [IU [12 and [131 From the theoretical point of view G(r, s) is the most fundamental 
property from which all others can be derived. It is presented last since it is computationally 
more demanding — at least if high accuracy is needed — and coefficients such as b e may 
be best determined directly from the moments. The normalized histograms G(r, s) are 
computed by counting the number of segment vectors between r — dr/2 and r + dr/2 with 
dr being the width of the bin and one divides then by the spherical bin volume. Since 
the BFM model is a lattice model, this volume is not 4nr 2 dr but given by the number of 
lattice sites the segment vector can actually point to for being allocated to the bin. Incorrect 
histograms are obtained for small r if this is not taken into account. (Averages are taken 
over all segments and chains, just as before.) Clearly, non-universal physics must show up 
for small vector length r and small curvilinear distance s and we concentrate therefore on 
values r ^> a and s > 31. 

When plotted in linear coordinates as in Fig. [TU G(r, s) compares roughly with the 
Gaussian prediction G (r, s) given by Eq. ([2]), but presents a distinct depletion for small 
segment sizes with n = rjb c ^fs <C 1 and an enhanced regime for n ~ 1. A second depletion 
region for large n 3> 1 — expected from the finite extensibility of the segments — can be best 
seen in the log- log representation of the data (not shown). To analyse the data it is better 
to consider instead of G(r, s) the relative deviation SG(r, s)/Go(r, s) = G(r, s)/Go(r, s) — 1 
which should further be divided by the strength of the segmental correlation hole, c e /\/s. 
As presented in Fig. [12] this yields a direct test of the key relation Eq. ([3]) announced in the 
Introduction. The figure demonstrates nicely the scaling of the data for all s and for both 
models. It shows further a good collapse of the data close to the universal function f(n) 
predicted by theory (bold line). Note that the depletion scales as 1/n for small segment 
sizes (dashed line). The agreement of simulation and theory is by all standards remarkable. 
(Obviously, error bars increase strongly for n ^> 1 where Go(r,s) decreases strongly. The 
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regime for very large n where the finite extensibility of segments matters has been omitted 
for clarity.) We emphasize that this scaling plot depends very strongly on the value b e which 
is used to calculate the Gaussian reference distribution. 

If a precise value is not available we recommend to use instead the scaling variable 
m = r/R e (s) for the horizontal axis, i.e. to replace the scale b c ^fs estimated from the 
behaviour of asymptotically long chains by the measured (mean-squared) segment size 
for the given s. The Gaussian reference distribution is then accordingly Go(m, R e (s)) = 
(3/2nR e (s) 2 ) exp(— |m 2 ). The corresponding scaling plot is given in Fig. [13j It is simi- 
lar and of comparable quality as the previous plot. Changing the scaling variable from 
n = r/b e y/s to m — r/R c (s) ~ (r/6 c i/i)(l + c e /2^fs) changes somewhat the universal func- 
tion. Expanding the previous result, Eq. ([3]), this adds even powers of m to the function 
f{n) given in Eq. ([3]) 

/(») => /H = (-^/f - i + 9m 4- yf m 2 - \m* Y (29) 

That the two additional terms in the function are correct can be seen by computing the 
second moment 4tt J drr A SG(r, s) which must vanish by construction. The rescaled relative 
deviation is somewhat broader than in the previous plot due to the additional term scaling 
as m 2 . As already stressed this scaling does not rely on the effective bond length b c and 
is therefore more robust. It has the nice feature that it underlines that there is only one 
characteristic length scale relevant for the swelling induced by the segmental correlation 
hole, the typical size of the chain segment itself. 

V. CONCLUSION 

Issues covered and central theoretical claims. We have revisited Flory's famous ideality 
hypothesis for long polymers in the melt by analyzing both analytically and numerically 
the segmental size distribution G(r, s) and its moments for chain segments of curvilinear 
length s. We have first identified the general mechanism that gives rise to deviations from 
ideal chain behavior in dense polymer solutions and melts (Sec. HI])- This mechanism rests 
upon the interplay of chain connectivity and the incompressibility of the system which 
generates an effective repulsion between chain segments (Fig. [2]). This repulsion scales like 
u(s) ~ c c / yfs where the "swelling coefficient" c e ~ l/blp se ^ s t ne strength of the interaction. 

23 



It is strong for small segment length s, but becomes weak for s — > N in the large- N limit. 
The overall size of a long chain thus remains almost 'ideal', whereas subchains are swollen as 
described by Eq. §5§. Most notably, the relative deviation SG(r, s)/G (r, s) of the segmental 
size distribution from Gaussianity should be proportional to u(s). As a function of segment 
size r, the repulsion manifests itself by a strong 1 /r-depletion at short distances r <C b e ^fs 
and a subsequent shift of the histogram to larger distances (Eq. (j3j)). 

Summary of computational results. Using Monte Carlo and molecular dynamics simu- 
lation of two coarse-grained polymer models we have verified numerically the theoretical 
predictions for long and flexible polymers in the bulk. We have explicitly checked (e.g., 
Figs. EE fTB"j) that the relative deviations from Flory's hypothesis scale indeed as 1/ y/s. 
Especially, the measurement of the bond-bond correlation function P(s), being the second 
derivative of the second moment of G(r, s) with respect of s, allows a very precise verifi- 
cation (Fig. |Sj) and shows that higher order corrections beyond the first-order perturbation 
approximation must be small. The most central and highly non-trivial numerical verification 
concerns the data collapse presented in Figs. Q21 and [13] for the segmental size distribution 
of both computational models. All other statements made in this paper can be derived and 
understood from this key finding. It shows especially that the swelling coefficient c e must 
be close to the predicted value, Eq. (j3J). 



It is well known 



lOj that the effective bond length is difficult to predict at low com- 



pressibility and no attempt has been done to do so in this paper. We show instead how the 
systematic swelling of chain segments - once understood - may be used to extrapolate for 
the effective bond length of asymptotically long chains. Figs. [5] and [6] indicate how this may 
be done using Eq. (jSJ). The high precision of our data is demonstrated in Fig. [121 by the 
successful scaling of the segmental size distribution. 

For several moments (r 2p ) we have also fitted empirical swelling coefficients c p using 
Eq. (jSJ). In contrast to the effective bond length b e these coefficients are rather well pre- 
dicted by one- loop perturbation theory if the bond length b of the reference Hamiltonian is 
renormalized to the effective bond length b e , as we have conjectured in Sec. Ill Bt Since the 
empirical swelling coefficients, c p ~ c e (b/b e ) 2p ~ 3 , would otherwise strongly depend on the 
moment taken, as shown in Eq. (fT2|) . our numerical data (Tab. [XT]) clearly imply b/b c m 1. 
Minor deviations found for the BSM samples may be attributed to the fact that monomers 
along the BSM chains do strongly overlap — an effect not taken into account by the theory. 
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To clarify ultimately this issue we are currently performing a numerical study where we 
systematically vary both the compressibility and the bond length of the BSM. 

General background and outlook. The most striking result presented in this work con- 
cerns the power law decay found for the bond-bond correlation function, P(s) oc 1/s 3 / 2 
(Fig. [S]). This result suggests an analogy with the well-known long-range velocity correla- 



tions found in dense fluids by Alder and Wainwright nearly fourty years ago [41], |42j. In 
both cases, the ideal uncorrelated object is a random walker which is weakly perturbed (for 
d > 2) by the self-interactions generated by global constraints. Although these constraints 
are different (momentum conservation for the fluid, incompressibility for polymer melts) the 
weight with which these constraints increase the stiffness of the random walker is always 
proportional to the return probability. It can be shown that the correspondence of both 
problems is mathematically rigorous if the fluid dynamics is described on the level of the 
linearized Navier-Stokes equations [43| . 

We point out that the physical mechanism which has been sketched above is rather 
general and should not be altered by details such as a finite persistence length — at least 
not as long as nematic ordering remains negligible and the polymer chains are sufficiently 
long. (Similarly, velocity correlations in dense liquids must show an analytical decay for 
sufficiently large times irrespective of the particle mass and the local static structure of 
the solution.) While this paper focused exclusively on scales beyond the correlation length 
of the density fluctuations, i.e. g( < 1 or s/g > 1, where the polymer solution appears 
incompressible, effects of finite density and compressibility can be readily described within 
the same theoretical framework and will be presented elsewhere 43] . To test our predictions, 
flexible chains should be studied preferentially, since the chain length required for a clear- 
cut description increases strongly with persistence length. This is in fact confirmed by 
preliminary and on-going simulations using the BSM algorithm. 

In this work we have only discussed properties in real space as a function of the curvilin- 
ear distance s. These quantities are straightforward to compute in a computer simulation 
but are barely experimentally relevant. The non-Gaussian deviations induced by the seg- 
mental correlation hole arise, however, also for an experimentally accessible property, the 
intramolecular form factor (single chain scattering function) F(q). As explained at the 
end of the Appendix, the form factor can be readily obtained by integrating the Fourier 
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transformed segmental size distribution given in Eq. ([3]). This yields 



q 2 F(q) 



12 

62 



1 - 



3 b e q 



(30) 



in agreement with the result obtained in Refs. [l4j, ll9| by direct calculation of the form factor 
for very long equilibrium polymers. As a consequence of this, the Kratky plot (q 2 F(q) vs. 
wave vector q) should not exhibit the plateau expected for Gaussian chains in the scale-free 
regime, but rather noticeable non-monotonic deviations. See Fig. 3 of [3]. This result 
suggests to revisit experimentally this old pivotal problem of polymer science. 

Our work is part of a broader attempt to describe systematically the effects of correlated 
density fluctuations in dense polymer systems, both for static [l2j, [l3j, |4J, |45j and dynam- 



ical 



29, 



35. 



46j| properties. An important unresolved question is for instance whether the 



predicted long-range repulsive forces of van der Waals type ( "Anti-Casimir effect" 
are observable, for instance in the oscillatory decay of the standard density pair-correlation 
function of dense polymer solutions. Since the results presented here challenge an important 
concept of polymer physics, they should hopefully be useful for a broad range of theoreti- 
cal approaches which commonly assume the validity of the Gaussian chain model down to 
molecular scales 47, Q, Q]- This study shows that a polymer in dense solutions should not 



be viewed as one soft sphere (or ellipsoid) [50|, l5ll, 152J , but as a hierarchy of nested segmental 



correlation holes of all sizes aligned and correlated along the chain backbone (Fig. El (b)). We 
note that simi 
polymers 



ar deviations from Flory's hypothesis have been reported recently for linear 
47 ] and polymer gels and networks 0, Q. 
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The repulsive interactions 



should also influence the polymer dynamics, since strong deviations from Gaussianity are 
expected on the scale where entanglements become important, hence, quantitative predic- 
tions for the entanglement length N e have to be regarded with more care. The demonstrated 
swelling of chains should be included in the popular primitive path analysis for obtaining N e 



551 ]. especially if 'short' chains (N < 500) are considered. The effect could be responsible 



p,q 



N 



(1 

1 |s, 



561 ] as may be seen by considering the cor- 
; (X p • X 3 ) of the Rouse modes X p = i J^ =0 dnr n cos(npir / N) where 
57|. Using (r n - r r , 



for observed deviations from Rouse behavior \2i 
relation function C n 

) 2 = rf l + rfn — 2r n ■ r m for the segment size, this 



correlation function can be readily expressed as an integral over the second moment of the 
segmental size distribution 

1 



a 



pq 



2N 2 



N 



N 



dn I dm((i 



J 2 ) cos(np7r / TV) cos(mp7r / TV) 



(31) 
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which can be solved using our result Eq. (j3j). This implies for instance for p = q that 



C pp 



Nbl 



1 



7T 



(32) 



The bracket entails an important correction with respect to the classical description given by 
the prefactor jf|. We are currently working out how static corrections, such as those for C pp , 
may influence the dynamics for polymer chains without topological constraints. (This may 
be realized, e.g., within the BFM algorithm by using the L26 moves described in Sec. IIII Al ) 
Moreover, for thin polymer films of width H the repulsive interactions are known to be 
stronger than in the bulk 12|. This provides a mechanism to rationalize the trend towards 



swelling observed experimentally 58] and confirmed computationally 21] : 



1 - 



bis 



hg(s)/H. 



(33) 



(Prefactors omitted for clarity.) Here R x (s) and b x denote the components of the segment 
size and the effective bond length parallel to the film. It also explains the (at first sight 
surprising) systematic increase of the polymer dynamics with decreasing film thickness 22]. 
Specifically, the parallel component of the monomer mean-squared displacement g x (t) is 
expected to scale as g x (t) ~ Rl(s(t)) oc t 1//4 (l + \og{t)/H) for long reptating chains where 
s{t) oc t 1//4 6]. (The corresponding effect for the three-dimensional bulk should be small, 
however.) For the same reason (flexible) polymer chains close to container walls must be 
more swollen and, hence, faster on intermediate time scales than their peers in the bulk. 
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APPENDIX A: MOMENTS OF THE SEGMENTAL SIZE DISTRIBUTION AND 
THEIR GENERATING FUNCTION 



Higher moments of the segmental size distribution G(r, s) can be systematically obtained 
from its Fourier transformation 



G(q,s) = J d 3 rG(r,s) exp(iq-r) 



which is in this context sometimes called the "generating function" 59]. For an ideal 
Gaussian chain, the generating function is then Go(q,s) = exp(— sq 2 a 2 ) where we have 
used a 2 = b 2 /6 instead of the bond length b 2 to simplify the notation. Moments of the 
size distribution are given by proper derivatives of G(q, s) taken at q = 0. For example, 
(r 2p ) = (— l) p A p G(q, s) q= o (with A being the Laplace operator with respect to the wave vec- 
tor q). A moment of order 2p is, hence, linked to only one coefficient A 2p in the systematic 
expansion, G(q, s) = J2 p=0 A 2p q 2p , of G(q, s) around q — 0. For our example this implies 

(r 2p ) = (-l) p (2p+l)\A 2p (Al) 

in general and more specifically for a Gaussian distribution (r 2p ) = ( - 2p p \ 1 ^' s p a 2p . The non- 
Gaussian parameters read, hence, 



6V A 2p 

Ms) -l {2p + 1} , (r2)p - 1 p. A ,, (A2) 



which implies (by construction) a p = for a Gaussian distribution. As various moments of 
the same global order 2p are linked to the same A 2p they differ by a multiplicative constant 
independent of the details of the (isotropic) distribution G(q, s). For example, (r 2 ) = 6|A 2 |, 
(r 4 ) = 120A 4 , (x 2 ) = (y 2 ) = 2\A 2 \, (x 2 y 2 ) = 8A 4 with x and y denoting the spatial 
components of the segment vector r. Using Eq. (1A2I) for p = 2 it follows that 

^W-i-|gy = i-2| = M S ), (A3) 

i.e. the properties a 2 (s) and K xy (s) discussed in Figs. 151 and [TU1 must be identical in general 
provided that G(q, s) is isotropic and can be expanded in q 2 . 

We turn now to specific properties of G(q, s) computed for formally infinite polymer chains 
in the melt. In practice, these results are also relevant for small segments in large chains, 
N ^> s ^> 1, and, especially, for segments located far from the chain ends. These chains are 
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nearly Gaussian and the generating function can be written as G(q, s) = Go(q, s) + SG(q, s) 
where SG(q, s) = —(UG)o + (U)q(G)o is a small perturbation under the effective interaction 
potential v(q) given by Eq. Q. To compute the different integrals it is more convenient to 
work in Fourier-Laplace space (q, t) with t being the Laplace variable conjugate to s: 



5G(q,t) 



ds SG(q, s)e 



-St 



As illustrated in Fig. [HI there a three contributions to this perturbation: one due to in- 
teractions between two monomers inside the segment (left panel), one due to interactions 
between an internal monomer and an external one (middle panel) and one due to interac- 
tions between two external monomers located on opposite sides (right panel). In analogy to 



the derivation of the form factor described in Ref. 

1 v 



14j this yields: 



5G(q,t) 



+ 



(q 2 a 2 + t) 2 Aira 3 

1 v 
{q 2 a 2 + 1) 2 Anqa 2 ^ 2 

I 2v 
q 2 a 2 + 1 Anqa 4 



q l a 2 + 1 



Vi) 



Arctan 


qa 


a/i + Vt_ 





qa 



Arctan 


qa 


a/i + y/t. 





a/i + ^q 2 a 2 + t / 
qa 



Airqa t 



Arctan 



qa 

ivi. 



qa 



(A4) 



y / q 2 a 2 + t / 

The graph given in the left panel of Fig. [14] corresponds to the first two lines, the middle 
panel to the third line and the right panel to the last one. Seeking for the moments we 
expand 5G(q, t) around q = 0. Having in mind chain strands counting many monomers 
(s ^> 1), we need only to retain the most singular terms for t — > 0. Defining the two 
dimensionless constants d = v^/3ira 4 = I2v£/7rb 4 and c = (37r 3 / 2 a 3 p) -1 = a/24/7t 3 /6 3 p this 
expansion can be written as 

ir(2) 



8G(q,t) 



+ 



1! t 2 

2 r(3) 

3 T(4) 
3!^ 



d a 2 q 2 + 

j 4 4 

d a q — 
d a 6 q 6 + 



1 T(3/2) 2 2 
1 16 1X5/2) 44 

2iy^ cag 

3! 35 f' 2 q + ' 



(A5) 



where we have used Euler's Gamma function T(a) 60|. The first leading term at each order 
in q 2 — being proportional to the coefficient d — ensures the renormalization of the effective 
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bond length. The next term scaling with the coefficient c corresponds to the leading finite 
strand size correction. Performing the inverse Laplace transformation T(a)/t a — > s"" 1 and 
adding the Gaussian reference distribution G (q,s) this yields the /^-coefficients for the 
expansion of G(q, s) around q = 0: 



.4 



o 



A 2 = -a 2 s il + d °-= 



A, = 


1 4 2 

-as 


A 6 = 





A, = 



16 c 
216 c 



35 y/s 



(A6) 



More generally, one finds 



^ = ^<^l 1+ ^-^n%75j (A7) 

From this result and using Eq. flAll) one immediately verifies that the moments of the 
distribution are given by the Eqs. (jTTj) and (|T2|) . Using Eq. (1A2|) one justifies similarly 
Eq. ( 1261) for the non- Gaussian parameter a p . 

These moments completely determine the segmental distribution G(r, s) which is indi- 
cated in Eq. ffT3"j) . While at least in principle this may be done directly by inverse Fourier- 
Laplace transformation of the correction 5G(q, t) to the generating function it is helpful to 
simplify further Eq. (1A4I) . We observe first that 5G(q,t) does diverge for strictly incom- 
pressible systems (v — > oo) and one must keep v finite in the effective potential whenever 
necessary to ensure convergence (actually everywhere but in the diagram corresponding to 
the interaction between two external monomers). Since we are not interested in the wave 
vectors larger than l/£ we expand 5G(q,t) for £ — > which leads to the much simpler 
expression 

The first term diverges as -^/w for diverging v. It renormalizes the effective bond length in 
the zero order term which is indicated in the first line of Eq. (|T3|) . The next two terms scale 
both as v°. Subsequent terms must all vanish for diverging v and can be discarded. It is 
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then easy to perform an inverse Fourier-Laplace transformation of the two relevant v terms. 
This yields 

= g '<^¥(-§ + t-t) < A9 > 

with x = rja^fs = \^6n. This is consistent with the expression given in the second line of 
Eq. m. 



We note finally that the intramolecular form factor F(q) = -k J2 n m=i ( ex P(^ ' ( r n ~ r m)) 
of asymptotically long chains can be readily obtained from Eq. ( 1A8I) . Observing that 
(exp(iq • (r n — r m )) = J <i 3 rexp(zq • r)G(r, s) = G(q, s) one finds 

5F(q) = 2 J ds 5G(q, s) = 2 5G(q, t = 0) = ~ 2 ^^^ ; ( A1 °) 

where we used the third term of Eq. (1A8|) in the last step. The first term in Eq. (1A8j) is 
discarded as before, since it renormalizes the effective bond length in the reference form 
factor: Fo(q) = 12/b 2 q 2 =^ 12/blq 2 . It follows, hence, that within first-order perturbation 
theory 

F(q) = F (q) + 6F(q) « F (q) (l - 3 -^j (All) 

as indicated by Eq. fl30|) in the Conclusion. This is equivalent to the result 1/F(q) — 1/F n (q) w 
q 3 /32p discussed in Refs. [14|, [l9| for polymer melts and anticipated by Schafer |11| by 



renormalization group calculations of semidilute solutions. 
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N 






R e 


Re 


b e (N) 


\/6bJN) 


16 


2 16 


1214 


11.7 


4.8 


2.998 


2.939 


32 


2 15 


3485 


17.1 


7.0 


3.066 


3.030 


64 


2 14 


1.1 • 10 4 


24.8 


10.1 


3.116 


3.094 


128 


8192 


3.3 • 10 4 


35.6 


14.5 


3.153 


3.139 


256 


4096 


1.0- 10 5 


50.8 


20.7 


3.179 


3.171 


512 


2048 


3.2 • 10 5 


72.2 


29.5 


3.200 


3.193 


1024 


1024 


1.0- 10 6 


103 


42.0 


3.216 


3.212 


2048 


512 


3.2 • 10 6 


146 


59.5 


3.227 


3.223 


4096 


256 


9.7- 10 6 


207 


85.0 


3.235 


3.253 


8192 


128 


2.9- 10 7 


294 


120 


3.249 


3.248 



TABLE I: Various static properties of dense BFM melts of number density p = 0.5/8: the chain 
length N, the number of chains n c h per box, the relaxation time r e characterized by the diffusion 
of the monomers over the end-to-end distance and corresponding to the circles indicated in Fig. [3l 
the root-mean-squared chain end-to-end distance R e and the radius of gyration R g of the total 
chain (s = N — 1). The last two columns give estimates for the effective bond length from the 
end-to-end distance, b e (N) = R e /(N — l) 1 / 2 , and the radius of gyration, b g (N) = R g /\fN. The 
dashed line in Fig. 2] indicates b e (N) 2 . Apparently, both estimates increase monotonicly with 
reaching b e (N) \/6b g (N) m 3.2 for the largest chains available. Note that \/6b g (N) < b e (N) for 
smaller N. 
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Property 


BFM 


BSM 


Length unit 


lattice constant 


bead diameter 


Temperature k^T 


1 


1 


Number density p 


0.5/8 


0.84 


Linear box size L 


256 


< 62 


Number of monomers n mon 


1048576 


< 196608 


Largest chain length N 


8192 


1024 


Mean bond length (|l n |) 


2.604 


0.97 


I = (12)V2 


2.636 


0.97 


Effective bond length b c 


3.244 


1.34 


Pb 3 e 


2.13 


2.02 


Coo = {b c /l) 2 


1.52 


1.91 


l P = l(Coo + l)/2 


3.32 


1.41 


c e = ^24/ir 3 /pb 3 e 


0.41 


0.44 


ci/c 


1.0 


1.2 


C2/C e 


1.0 


1.1 


C3/c e 


1.0 


1.0 


C4/c e 


1.1 


1.2 




1.1 


0.9 


cp = ci(6 e /Z) 2 /8 


0.078 


0.124 


Dimensionless compressibility g 


0.245 


0.08 


Compression modulus v = 1/gp 


66.7 


14.9 


G z = y/vp/blp 


0.96 


1.8 



TABLE II: Comparison of some static properties of dense BFM and BSM melts. The first six rows 
indicate conventions and operational parameters. The effective bond length b c and the swelling 
coefficients c p (defined in Eq. (JSJ) are determined from the first five even moments of the segmental 
size distribution. The dimensionless compressibility g = S(q — > 0)/p has been obtained from the 
total static structure factor S(q) = Ylkl=i (exp(iq • (r^ — r;))) in the zero wave vector limit 
as shown at the end of Ref. [141 ] . The values indicated correspond to the asymptotic long chain 
behavior. Properties of very small chains deviate slightly. 
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FIG. 1: (Color online) Sketch of a polymer chain of length N in a dense melt in d = 3 dimensions. 
As notations we use rj for the position vector of a monomer i, lj = — rj for its bond vector, 
r = r m — r n for the end-to-end vector of the chain segment between the monomers n and m = n + s 
and r = ||r|| for its length. Segment properties, such as the 2p-th moments (r 2p ), are averaged over 
all possible pairs of monomers (n, m) of a chain and over all chains. The second moment (p = 1) 
is denoted R e (s) = (r 2 ) 1 ^ 2 , the total chain end-to-end distance is R e (s = N — 1). The dashed 
lines show the relevant graphs of the analytical perturbation calculation outlined in Sec. Ill Bi 
The numerical factors indicate for infinite chains (without chain end effects) the relative weights 
contributing to the l/y^-swelling of R e (s) indicated in Eq. (fTO]) . 
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FIG. 2: (Color online) Role of incompressibility and chain connectivity in dense polymer solutions 
and melts, (a) Sketch of the segmental correlation hole of a marked chain segment of curvilin- 
ear length s. Density fluctuations of chain segments must be correlated, since the total density 
fluctuations (dashed line) are small. Consequently, a second chain segment feels an entropic re- 
pulsion when both correlation holes start to overlap, (b) Self-similar pattern of nested segmental 
correlation holes of decreasing strength u(s) ~ s/pR(s) 3 c e /^/s aligned along the backbone of 
a reference chain. The large dashed circle represents the classical correlation hole of the total 
chain (s ~ N) 5|. This is the input of recent approaches to model polymer chains as soft spheres 



50, 



521 ] . We argue that incompressibility on all scales and chain connectivity leads to a short 



distance repulsion of the segmental correlation holes, which increases with decreasing s. 



35 




4 



FIG. 3: (Color online) Diffusion time r e over the (root-mean-squared) chain end-to-end distance 
R C (N — 1) as a function of chain length N for different versions of the Bond Fluctuation Model 
(BFM). All data indicated are for the high number density (p = 0.5/8) corresponding to a polymer 
melt with half the lattice sites being occupied. We have obtained r e = R^(N — 1)/6D S from the self- 
diffusion coefficient D s measured from the free diffusion limit of the mean-squared displacement 
of all monomers (5r(t) 2 ^ = QD s t. Data from the classical BFM with topology conserving local 
Monte Carlo (MC) moves in 6 spatial directions (L06) are represented as stars. All other data 
sets use topology violating local MC moves in 26 lattice directions (L26). If only local moves are 
used, L26-dynamics is even at relatively short times perfectly Rouse like which allows the accurate 
determination of D s although the monomers possibly have not yet moved over R e (N — 1) for the 
largest chain lengths considered. Additional slithering snake (SS) moves increase the efficiency of 
the algorithm by approximately an order of magnitude (squares, bold line). The power law exponent 
is changed from 2 to an empirical 1.62 (dashed line) if in addition we perform double-bridging (DB) 
moves. 
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FIG. 4: (Color online) Mean-squared segment size R e (s) 2 /s vs. curvilinear distance s. We present 
BFM data for different chain length N at number density p = 0.5/8. The averages are taken over 
all possible monomer pairs (n,m = n + s). The statistics deteriorates, hence, for large s. Log-linear 
coordinates are used to emphasize the power law swelling over several orders of magnitude of s. 
The data approach the asymptotic limit (horizontal line) from below, i.e. the chains are swollen. 
This behavior is well fitted by Eq. ((SJ) for 1 <C s -C N (bold line). Non-monotonic behavior is found 
for s — > N, especially for small N. The dashed line indicates the measured total chain end-to-end 
distances, b e (N) 2 = R e (N - 1) 2 /{N - 1) from Tab. U showing even more pronounced deviations 
from the asymptotic limit. The dash-dotted line compares this data with Eq. (|19|) . 
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FIG. 5: (Color online) Replot of the mean-squared segment size as y = K%(s) = 1 — R c (s) 2 /b%s 
vs. x = ci/i/s, as suggested by Eq. ([5]), for different trial effective bond lengths b c as indicated. 
Only BFM chains of length N = 2048 are considered for clarity. This procedure is very sensitive 
to the value chosen and allows for a precise determination. It assumes, however, that higher order 
terms in the expansion of K±(s) may be neglected. The value b e is confirmed from a similar test 
for higher moments (Fig. [6]) . 
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FIG. 6: (Color online) Critial test of Eq. © where the rescaled moments y = K p (s) of the segment 
size distribution (defined in Eq. (pQ)) are plotted vs. x = 2(2p+i)! ^/f- We consider the first five even 
moments (p = 1, . . . , 5) for the BFM with N = 2048 and the BSM with N = 1024. Also indicated 
is the rescaled radius of gyration, y = 5/8 (1 — 6Rg(s)/b^,(s + 1)), as a function of x = c\ /\/s + 1 
(filled circles). The BSM data has been shifted upwards for clarity. Without this shift a perfect 
data collapse is found for both models and all moments. Keeping the same effective bond length 
b e for all moments of each model we fit for the swelling coefficients c p by rescaling the horizontal 
axis. We find b e ~ 3.244 for the BFM and 1.34 for the BSM. If b e is chosen correctly, all data 
sets extrapolate linearly to zero for large s (x — ► 0). The swelling coefficients found are close the 
theoretical prediction c e , as indicated in Tab. HH The plot demonstrates that the non-Gaussian 
deviations scale as the segmental correlation hole, u(s) ~ c e / 'y/s and this for all moments as long 
as x -C 1. The saturation at large x is due to the finite extensibility of short chain segments. Since 
this effect becomes more marked for larger moments, the fit of b c is best performed for p = 1. 
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FIG. 7: (Color online) Plot of K\(s) as a function of u{s)c\/c e ~ l/yfs using the measured 
u(s) = y / 24/7r 3 s/pi? c (s) 3 . For A = 2 (corresponding to two segments being connected) BFM and 
BSM data are compared. Several A values are given for N = 2048 BFM chains. For chain segments 
with 1 -C s <C N all data sets collapse on the bisection line confirming the so-called "recursion 



relation" K\ ~ u proposed by Semenov and Johner 12]. The statistics becomes insufficient for 
large s (left bottom corner). Systematic deviations arise for s — > N due to additional finite- N 
effects. 
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FIG. 8: (Color online) The bond-bond correlation function P(s)/cp as a function of the curvilinear 
distance s. Various chain lengths are given for BFM. Provided that 1 <C s <S N, all data sets 
collapse on the power law slope with exponent uj = 3/2 (bold line) as predicted by Eq. (|23p . The 
dash-dotted curve P(s) ~ exp(— s/1.5) shows that exponential behavior is only compatible with 
very small chain lengths. The dashed lines correspond to the theoretical prediction, Eq. (|24p . for 
short chains with N = 16, 32, 64 and 128 (from left to right). 
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FIG. 9: (Color online) Non-Gaussian parameter a p (s) computed for the end-to-end distance of 
chain segments as a function of Ce/y/s. Perfect data collapse for all chain lengths and both sim- 
ulation models is obtained for each p. A linear relationship over nearly two orders of magnitude 
is found as theoretically expected. Data for three moments (p = 2, 3, 4) are indicated showing 
a systematic increase of non-Gaussianity with p. The data curvature for small s becomes more 
pronounced for larger p. 
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FIG. 10: (Color online) Plot of K xy (s) = 1 — (x 2 y 2 ^ / (x 2 ) (y 2 ) averaged over all pairs of monomers 
(n,m = n + s) and three different direction pairs as a function of c c /yfs. As indicated by the sketch 
at the bottom of the figure, K xy {s) measures the correlation of the components of the segment 
vector r. All data points collapse and show again a linear relationship K xy ~ u{s). Different 
directions are therefore coupled! No curvature is observed over two orders of magnitude confirming 
that higher order perturbation corrections are negligible. Noise cannot be neglected for large 
s > 100 and finite segment-size effects are visible for s ~ 1. 
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FIG. 11: (Color online) Segment size distribution y = G(r, s)(b e s 1 / 2 ) 3 vs. n = r/beS 1 ! 2 for several 
s as indicated in the figure. Only data for BFM with N = 2048 and BSM with N = 1024 are 
presented. (A similar plot can be achieved by renormalizing the axes using R e {s) instead of freS 1 / 2 ). 
The bold line denotes the Gaussian behaviour y = (3/2-/r) 3 / 2 exp(— 3n 2 /2). One sees that compared 
to this reference the measured distributions are depleted for small n <C 1 (where the data does not 
scale) and enhanced for n ~ 1. 
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Enhancement 
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FIG. 12: (Color online) Deviation 5G(r,s) = G(r,s) — Go(r,s) of the measured segmental 
size distribution from the Gaussian behavior Go(r,s) expected from Flory's hypothesis for sev- 
eral s and both models as indicated in the figure. As suggested by Eq. (|3|), we have plotted 
y = (5G(r,s)/Go(r,s))/(c e /y/s) as a function of n = r/b e yfs. The Gaussian reference distribution 
has been computed according to Eq. ([2]) for the measured effective bond length b e . A close to 
perfect data collapse is found for both models. This shows that the deviation scales linearly with 
u(s) « Ce/y^, as expected. The bold line indicates the universal function of f(n) predicted by 
Eq. ©. 
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m=r/R (s) 

FIG. 13: (Color online) Replot of the relative deviation of the measured segment size distribution, 
y = (SG(r, s)/Go(r, s)))/(c e /v / s), as a function of m = r/R e (s). The figure highlights that the 
measured segment size is the only length scale relevant for describing the deviation from Flory's 
hypothesis. The same data sets and symbols are used as in the previous Fig. [T2J 
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FIG. 14: Interaction diagrams used in reciprocal space for the calculation of 5G(q,t) in the scale 
free limit. There exist three nonzero contributions to first-order perturbation, the first involving 
two points inside the segment (first two lines of Eq. (|A4|) ). the second one point inside and one 
outside the segment (third line of Eq. (|A4j) ) and the third one point on either side of the segment 
(last line of Eq. (|A4j) ). Momentum q flows from one correlated point to the other. Integrals are 
performed over the momentum k. Dotted lines denote the effective interactions v(k) given by 
Eq. ©, bold lines the propagators which carry each a momentum q or q — k as indicated. 
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